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Abstract 


The voltammetry for the reduction of oxygen at a microdisc electrode is reported in two room 
temperature ionic liquids: 1-butyl-1-methylpyrrolidinium &*s(trifluoromethylsulfonyl)imide 
([Bmpyrr] [NTf 2 ]) and trihcxyltetradecylphosphonium 6z.s(trifluoromethylsulfonyl)imide 
([Pi 4 , 6 , 6 , 6 ][NTf 2 ]) at 298 K. Simulated voltammograms using Butler-Volmcr theory and symmet¬ 
ric Marcus-Hush (SMH) theory were compared with experimental data. Butler-Volmer theory 
consistently provided experimental parameters with a higher level of certainty than SMH the¬ 
ory. A value of solvent reorganisation energy for oxygen reduction in ionic liquids was inferred 
for the first time as 0.4-0.5 eV, which is attributable to inner-sphere reorganisation with a neg¬ 
ligible contribution from solvent reorganisation. The developed Butler-Volmcr and symmetric 
Marcus-Hush programs are also used to theoretically study the possibility of kinetically lim¬ 
ited steady state currents, and to establish an approximate equivalence relationship between 
microdisc electrodes and spherical electrodes resting on a surface for steady state voltammetry 
for both Butler-Volmcr and symmetric Marcus-Hush theory. 

1 Introduction 

Current understanding of interfacial electron transfer (ET) kinetics is dominated by two theo¬ 
ries, Butler-Volmer-- and Marcus-Hush,-^ which both describe a one-electron transfer: 

A + e-^B (1) 
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The first, Butler-Volmer theory (BV, Equations [2] and [3]), is well-reported and reliable for 
parameterisation,- but is phenomenological, and therefore is unable to provide any detailed 
physical insights of a predictive nature. BV theory relates the electrode kinetics to the hetero¬ 
geneous rate constants ( k °), the formal potential (Ef), and the transfer coefficient, a, which 
empirically describes how product or reactant-like the transition state is. 6 


k 


OX 


hed = A: 0 exp 



= A;°exp 


(1 -«) 
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RT 
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( 3 ) 


An alternative theory, symmetric Marcus-Hush (SMH), offers greater microscopic physical 
insight into the system under study. MH theory relates the heterogeneous rate constant to the 
reaction Gibbs energy, A G^, and reorganisation energy, A: 


A G x 



k° 


A exp 


'-AG*' 

RT 


( 4 ) 

( 5 ) 


As is shown in Equation |6l the total reorganisation energy is the sum of contributions from 
inner sphere bond reorganisation energy (A;) and solvent reorganisation energy (A 0 ). These can 
be estimated separately, with the outer-sphere solvent reorganisation energy from the Born- 
Marcus solvation equation,- and inner-sphere contributions from computational calculations 
(for example, DFT eto^) of changes in bond lengths and angles of the molecule. 


A — Ai + A 0 


( 6 ) 
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These two models have been used to analyse a range of electrochemical processes by cyclic 
voltammetry in aqueous^ 1 ^ and non-aqueous-AR— solvents. A comprehensive review by Hen- 
stridge et al — concluded that SMH theory struggles to adequately describe the voltammetric 
response of solution-phase redox couples whose transfer coefficients deviate significantly from 
0.5. A fourth parameter, 7 , was introduced, which accounts for any differences between the 
inner-shell force constant of the reactant and the product, to create asymmetric Marcus-Hush 
theory (AMH). This refinement was explored in aqueousA^ and non aqueous^ - — solvents and 
surface bound systems^ by Laborda et al, and it was concluded that AMH was able to pararn- 
eterise the experimental voltammetry at least as well as BV theory, and more accurately than 
SMH theory. 

The electroreduction of oxygen has been extensively studied in a variety of solvents.— 
The product of the initial one-electron reduction, superoxide, is a highly nucleophilic ion that 
initiates further reactions with a proton source, if one is present in solution. However, in 
absence of a proton source, the superoxide ion is stable, and the reaction proceeds as a simple 
one electron transfer: 

O 2 + e ^ O* (7) 

ET kinetics of irreversible oxygen reduction has been studied in-depth in aqueous solvents^— 
and values for the total reorganisation energy for the reduction of oxygen to superoxide have 
been estimated to lie between 1.26-4.51 eV. 2 ' This has been attributed to rearrangement of the 
polar water molecules upon formation of the superoxide ions. 

Room temperature ionic liquids (RTILs) are molten salts with melting points below room 
temperature,— and normally consist of a bulky, asymmetric organic cation, and an inorganic 
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anion.— Ionic liquids are interesting solvents in which to examine the reversible one-electron 
reduction of oxygen, because the superoxide has been known to be stable, leading to quasi- 
reversible voltammetry in a range of ILs.— >24i2!h22i 

Lynden-Bcll has conducted a thorough study into the applicability of Marcus theory to 
RTIL solvents.-^ 3 ' 34 It was suggested via simulation that the outer sphere reorganisation energy 
for redox processes in RTILs will be similar to that for the same process in a non ionic polar 
solvent. In a polar solvent, the molecules will re-orient themselves in response to the addition 
of charge, and in an ionic liquid the ions will translate themselves. It was predicted that these 
two processes produce similar reorganisation energies. 

The mechanism and ET kinetics have been elucidated for oxygen reduction to hydrogen 
peroxide in protic ionic liquids,— and in hydrophobic RTILs in the presence of water.— As yet, 
the ET kinetics of the one electron reduction of oxygen to superoxide in dry ionic solvents has 
not been studied, and no values for reorganisation energy for this process have been established. 

This study is the first to examine SMH theory of one electron oxygen reduction in ILs at 
microdisc electrodes. The structures of the ions in the ILs used in this study are shown in Fig. 
HI The present work aims to develop the theory of SMH voltammetry at microdisc electrodes 
and experimentally determine a value for solvent reorganisation, A, for oxygen reduction in two 
ionic liquids. The present work also aims to compare the ability of SMH theory and BV theory 
to simulate experimental reduction of oxygen and quantify a goodness of fit for each. This 
will allow conclusions to be reached about the relative ability of each theory to parameterise 
experimental data. 
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2 Theory 


We consider a one electron reduction of an electroactive species: 

A + e“ ^ B (8) 

at a microdisc electrode, using both Butler-Volnrer theory and symmetric Marcus-Hush theory 
to describe the electron transfer kinetics. Cyclic voltammetry is simulated using both models, 
and both are subsequently compared to experimental data (see section E]). 

This study utilises microdisc electrodes, which offer several advantages over macroelectrodes. 
Firstly, as the capacitance of the double layer is proportional to the size of the electrode, a con¬ 
sequence of the reduced size of the electrode will be naturally lower capacitance. Secondly, the 
small currents passed reduce the ohmic drop and allow a two-electrode system to be employed.^ 
Lastly, the smaller area enables a very small volume of ionic liquid to be studied (in a T-cell 
system, as described previously—), which is desirable for logistical reasons (expense, ability to 
dry the aliquot of IL more quickly). 

Cylindrical coordinates are employed to model the two dimensional microdisc electrode and 
simulate cyclic voltammetry. Fig. [2] shows a schematic diagram of the microelectrode (a) 
and the two dimensional cylindrical space used in our model (b). Also defined are the spatial 
coordinates r and z, as well as the electrode radius r e . In cylindrical space, Fick’s second law 
for diffusional mass transport is: 


dci f <9 2 Cj d 2 Ci 1 f)cj \ 

dt 1 \ dr 2 dz 2 r dr J 


( 9 ) 
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All symbols are defined in Table [1] The mass transport is considered to take place under full 
support (as is usually assumed when ionic liquids are used as solvents), and therefore electrical 
migration can be neglected. Further, the absence of mechanical stirring and any significant 
temperature gradients eliminate the presence of forced and natural convection respectively, 
allowing a diffusion only model to be used. 

The electron transfer reaction is only considered to occur at the electrode surface (r < r e , 
z — 0). The rate of this reaction is dependent on the potential applied to the electrode, E (V). 
In cyclic voltammetry, the potential is given by the equation: 

E = \E S — E v — ut\ + E v (10) 

where E s and E v are respectively the initial applied potential and the vertex potential (V), u is 
the scan rate (V s” 1 ) and t is the time since the start of the experiment (s). The flux of species 
A across the electrode due to electrolysis is then described as: 

-Da^ 7 - = k red c A - /coxCb ( 11 ) 

where the rate constants k re d and k ox are functions of the applied potential E, and c° A and Cg are 
the concentrations, at the electrode surface, of species A and B respectively. The mathematical 
forms of k red and k ox depend on whether Butler-Volrner or Marcus-Hush kinetics are employed, 
as detailed below. 
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2.1 Butler-Volmer Kinetics 

Within the Butler-Volmer formalism of electron transfer kinetics, k re d and k ox are considered 
to be exponentially dependent on the applied potential:— 


k ie d = k°ex p \—ol6\ 

(12) 

k ox = fc°exp [(1 — a) 9\ 

(13) 


where k° is the standard heterogeneous electrochemical rate constant (ms 1 ), a is the transfer 
coefficient and 6 is a dimensionless potential, defined as: 

e= W {E ~ S ?) ( 14 ) 

where Ef is the formal potential of the A/B couple and F, R , and T have their usual meanings. 

The transfer coefficient a, as described in the introduction, has a value of between zero and 
one, and is considered to reflect the structure of the transition state of the A/B redox couple.- 


2.2 Marcus-Hush Kinetics 

In (symmetric) Marcus-Hush theory,—~ 46 k ie d and k ox are functions of the applied electrode 
potential and the reorganisation energy, A. Mathematically: 


krod k 


L = k ' 


0 Sred ($> 


Sred (0, A) 

0 S OX ( 9 , A) 
'S'ox (0, A) 


(15) 


(16) 




Sred/ox (#, A) is given by the following integral: 


‘S'red / ox (M) = 


exp 


-AG*, /RT 

red/ox/ 


1 + exp [=pcc] 


—da; 


(17) 


AGj a|/ox is the activation energy of reduction/oxidation, and is given by: 


agI 


d/ox A / 0 ~h X 

A 


RT 4 


(18) 


where A is the dimensionless reorganisation energy, given by: 


A = -^-A 

RT 


(19) 


x is an integration variable, defined as: 




( 20 ) 


This accounts for the continuum of energy levels in the electrode of energy e. 


2.3 Simulation Procedure 

To simulate cyclic voltammetry, the mass transport equation and electron transfer kinetic equa¬ 
tions are normalised to remove scaling factors from the simulation. Substituting dimensionless 
parameters as defined in Table [2] into the mass transport equation, we obtain: 


dC 1 = D ' ( d 2 G 1 OCA 

dr l \dR 2 + dZ 2 + RdRJ 


(21) 
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The boundary condition at the electrode surface (R < 1, Z — 0) becomes: 


dC A 

dZ 


K red C° A - K m C% 


L ox^B 


( 22 ) 


where 


Jlred = Jl°exp \—a9\ 

(23) 

K ox = K°ex p [(1 - a) 6} 

(24) 


for Butler-Volmer kinetics, and 


K m] = K 


K m = K 


p S'red (0, A) 
S red (0,A) 
o S ox ( 6 , A) 
S'ox (0, A) 


(25) 

(26) 


for Marcus Hush kinetics. In both cases the dimensionless standard heterogeneous electrochem¬ 


ical rate constant is defined as: 


K° 


k°r e 


(27) 


At the electrode surface, equal and opposite flux acts as the boundary condition for species 


B: 


dC B dC A 

B dZ dZ 


(28) 
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The outer boundaries of the simulation space are set at: 


-f^max 1 6 \/^max^niax 

(29) 

Zmax 6'\/-D max T m ax 

(30) 


where D' max and r max are, respectively, the largest dimensionless diffusion coefficient in the 
system and the total dimensionless time of the experiment. These boundaries have been shown 
to be sufficiently distant from the electrode to be well outside the depletion zone.— - — 

At these boundaries, as well as at the symmetry boundary at R — 0, all Z and at the 
insulating surface (R > 1, Z — 0), a zero flux condition is imposed: 


dC\ 

m 


o 


(31) 


where N is one of the spatial coordinates R or Z as appropriate. 

Initial conditions are C A = 1 and C-q = 0 throughout the simulation space. At r = 0 
the experiment starts and the electrode surface boundary conditions are imposed. The total 
dimensionless flux J at the electrode is calculated as: 

j= H^L d * (32) 

This dimensionless flux is transformed into a real current / by: 


/ = 2nFD A c* A r e J 


(33) 
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2.4 Homogeneous kinetics 


In more complex reaction mechanisms, the single electron transfer thus far simulated may be 
followed by homogeneous kinetics. For example, the reduction of oxygen in the ionic liquid 
[P 14 , 6 , 6 , 6 ] [NTf 2 ] (see Section 15.21) follows the following mechanism: 


or x 

x ^ Y + e" 


(34) 

(35) 

(36) 


where the homogeneous reaction involves proton abstraction by the superoxide ion from the 
phosphonium cation in the ionic liquid. This irreversible, pseudo first order reaction forms a 
product which can then be oxidised, at less negative potentials than the initial electron transfer 
to oxygen. The homogeneous step is considered to have a fast rate constant. 

This being the case, if only the forward, reductive, wave of a cyclic voltammogram is 
measured and simulated, species X and Y may be neglected and only the first 2 steps in the 
above mechanism considered. To incorporate the homogeneous step into the simulations, the 
dimensionless mass transport equation is amended to: 


dC B 

dr 



d 2 C B d 2 C B 1 dC B \ 
dR 2 + dZ 2 + R~dR ) 


— kc b 


(37) 


where K is a dimensionless homogeneous rate constant, defined as: 


K = 



(38) 
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and k is a first order homogeneous rate constant for the reaction of species B. 

2.5 Computational Methods 

The above model is discretised using the Crank-Nicolson method— and solved numerically over 
discrete spatial and temporal grids. The temporal grid is defined initially in terms of the applied 
dimensionless potential, 6. A parameter Ng is defined as the number of discrete time steps per 
unit 9. The value of 6 at each time step k is then defined: 

e k = e k _ x t ( 39 ) 

tv g 

where the minus sign is used for the reductive sweep, and the plus sign for the oxidative sweep. 
The dimensionless time, r, at each time step is then calculated as: 

T k = Tfc_ 1 + -J— (40) 

Ngcr 

A non regular, exponentially expanding spatial grid is used to ensure computational effi¬ 
ciency. In the Z direction, the first spatial point (at the clectrode/insulating surface) is defined 
as Z G = 0. The first step size is defined as A, with subsequent step sizes growing increasingly 
large until the simulation boundary is reached. Each Z direction spatial point is defined as: 

Zj = Zj _i + 7 (Zj _i — Zj_ 2 ) (41) 

The R direction spatial grid is defined in the same way, except that rather than expanding 
from a simulation boundary, the grid expands away from the edge of the electrode at R = 1 in 
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both directions until R — 0 and R = R, 


max 


are reached. Fig. [3] shows a schematic diagram of 


some of the spatial grid points used in simulations. 

In order to achieve convergence to within 0.5% of a fully converged result, the following 
grid parameters were used: Ng = 1000, A = 8 x 10 -5 and 7 = 1.25. The Alternating Direction 
Implicit (ADI) 51 method was employed to simulate the voltammetry, in conjuction with the 


Thomas Algorithm— to solve the mass transport equation and boundary conditions simultane¬ 


ously in matrix form. The model was coded in C++ with multithreaded parallel programming 
employed to increase efficiency. Typical running times were c.a. 10 minutes on a 2.26 GHz 
Intel(R) Xenon(R) CPU with 2.25 GB of RAM. 

2.6 Simulating voltammetry at impacting spherical nanoparticles 

Electrochemical reactions can occur at the surface of nanoparticles impacting upon a conduc¬ 
tive, but electrochemical inert surface. A schematic diagram of a spherical nanoparticle in 
contact with a conducting but otherwise inert surface is shown in Fig. [4] (a). Upon contact, the 
nanoparticle becomes electrochemically active, and a species in solution can be reduced/oxidised 
on its surface, if a suitable potential is applied to the inter surface. This is modelled in the 
same way as for a microdisc described above, but with changed boundary conditions.— - — A 
cylindrical coordinate system is still used, as defined in Fig. [4j Part (a) shows a schematic of 
the spherical electrode on the surface. Part (b) shows the 2 dimensional simulation space. 

The electrode surface boundary condition is now applied at over the line R 2 + (Z — l) 2 = 1 
where the flux normal to the electrode surface is given by: 



(42) 
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where N is a coordinate normal to the electrode surface, and k re and k ox are as defined above 
for either Butler-Volmcr or Marcus-Hush kinetics. Other boundary conditions were zero flux 
conditions at the R = 0,R = R max (= l + VAl^w), Z = 0 and Z = Z max (= 2 +vCD max r max ). 
The dimensionless flux is then defined as: 


J 


j-ir/2 

r dC A , dC A ] 

J—n/2 

dR + dZ Sm4, \ 


cos (f) d(p 


The current is then defined, as for a microdisc above, as: 


(43) 


/ = 2ttFD A c* A r e J (44) 

The temporal grid for the simulation of voltammetry at a spherical electrode is the same 
as defined above for a microdisc electrode. To define the spatial grid, the geometry of the 
electrode is taken into account. The electrode surface is divided into equal spaced points, 
with a grid vertex at each of these points. This defines the grid points in the region R < 1 and 
Z < 2, as shown schematically in Fig. |4j In the regions R > 1 and Z > 2, the grid expands 
exactly as outlined for the microdisc above. Convergence studies found = 200 sufficient to 
ensure convergence. Computational methods were identical to those employed for microdisc 
simulations. Simulation running times were c.a. 10 mins on a 2.26 GHz Intel(R) Xenon(R) 
CPU with 2.25 GB of RAM. 
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2.7 Chronoamperometry 


The Shoup and Szabo equation is used to calculate expected chronoamperometric responses at 
microdisc electrodes to within 0.6% error. 5 ' For a reduction: 


I = -4nFD A c A r e f (r) 


(45) 


/(r) = 0.7854 + 0.4432 t - 0 ' 5 + 0.2146exp [—0.3912r -0 ' 5 ] (46) 


r 



(47) 


2.8 Quantifying Error 

In order to examine the relative merits of Butler-Volmer theory and Marcus-Hush theory in 
simulating cyclic voltammetry, a quantitative measure of error is desirable. This study uses 
the Mean Scaled Absolute Deviation (MSAD) to quantify the error between experiment and 
theory. The MSAD value for a simulated fit is given by: 


MSAD 


N 


E 

N 


l exp 


^sir 


£ exp 


x 100% 


(48) 


where N is the number of experimental data points. Since it is unlikely that each experimental 
data point will have an exactly corresponding simulated data point, a linear interpolation is used 
between the two simulated data points surrounding the experimental data point in question. 
This interpolated simulated current is then used in the calculation of the MSAD value. 

Small values of the experimental current can lead to extremely large MSAD values despite 
the absolute difference between experiment and theory being small. For this reason, all exper- 


16 





imental data points with a value less than 5% of the maximum absolute value of the current 
are neglected. This includes the initial part of the reductive wave as well as a narrow part of 
the oxidative wave, where the current passes from negative to positive values. 

3 Theoretical results 

Previous work by Feldberg 45 and Henstridge et. al— showed that in the Marcus-Hush model 
of electron transfer kinetics, steady state currents can be kinetically rather than diffusionally 
limited for small values of the reorganisation energy, A. In Feldberg’s study, an approximate 
equation of Oldham and Zoski— describing steady state voltammetry at a microdisc was invoked 
and modified to use Marcus-Hush, rather than Butler-Volrner kinetics. Henstridge used numer¬ 
ical simulation to generate Marcus-Hush steady state responses at a conductive nanosphere 
resting on an electrochemically inert but conductive surface. In this study, we simulate steady 
state voltammetry numerically using the symmetric Marcus-Hush formalism of electron trans¬ 
fer, and compare the circumstances in which a kinetically limited current is observed at a 
microdisc to those at a nanosphere. We also extend previous work showing the equivalence 
of microdiscs and spheres-on-surfaces under the Butler-Volrner formalism (subject to a mass 
transport correction factor),— and show that the approximate equivalence is retained when 
Marcus-Hush kinetics are used. 

3.1 Kinetically limited currents at a microdisc electrode 

Steady state voltammetry at a microdisc electrode was simulated using the model outlined 
above. Various values of dimensionless reorganisation energy, A, and dimensionless heteroge- 
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neous rate constant, K°, were simulated, and the results are shown in Fig. [5j Part (a) shows 
results for K° = 5 x 10 -1 , (b) shows K° = 5 x 10~ 3 and (c) shows K° = 5 x 10 -5 . Also shown in 
each case is the steady state voltammetry simulated using Butler-Volrner kinetics with the same 
values of K°, and a = 0.5, and Marcus-Hush results are shown relative to this Butler-Volrner 
result. For all three values of K°, it is seen that large A values produce the same limiting cur¬ 
rent as the Butler-Volrner model. As A becomes smaller, the limiting current becomes less than 
predicted by the Butler-Volrner model as kinetic control takes over. This effect is more extreme 
for the smallest values of K°, where even a A value of 200 has not reached the Butler-Volrner 
limit at all points on the wave. For K° = 5 x 10 _1 , unfeasibly small A values would be required 
to show significant deviation form the Butler-Volrner result. This is in qualitative agreement 
with the previous studies of this nature. 45 ’ 55 

The definition of K 0 in our dimensionless model is: 


A' 0 


r e k° 


(49) 


Lowering K° will therefore mean that either k° or r e is lowered, or D a is raised. Lowering 
k° trivially means the electrochemical kinetics become more likely to be the limiting factor in 
the steady state current as opposed to diffusion. Either lower r e or increasing D a also both 
lower K° and lead to a more kinetically limited current. This is due to the fact that either 
of these changes will increase mass transport efficiency to the electrode surface, thus making 
electrochemical kinetics the more limiting factor. Therefore, in order to see a kinetically limited 
current, it is desirable to have a small k°, a small r e and a large D\. 
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3.2 Comparison of Marcus-Hush theory at microdisc and impacting 


nanoparticles 

Recent work by Molina et. al, 54 shows that steady state voltammetry at an isolated spherical 
electrodes alone in solution, a spherical electrode supported on a conductive but inert surface 
and a microdisc electrode (all shown schematically in Fig. [6]) may all be approximately described 
(under the Butler-Volmer formalism) by a single equation: 

I _ iCexp (-aO) 

him 1 + A" 0 exp (—a&) [1 + exp (6 1 )] 

where K° is a dimensionless heterogeneous rate constant scaled by a 
transport factor. These parameters are defined for isolated spheres 
surfaces (sos) and microdiscs (md) as: 


(50) 

geometry dependent mass 
(is) spheres supported on 


K° = 

IS 

k° r e 

Oa = A " 

(51) 

K° = 

SOS 

k°r, 1 

(52) 

Da ln2 ln2 

K° - 

n md — 

k°r e 7r AT° d 7T 

(53) 

Da 4 4 


Therefore, steady state voltammetry at a microdisc of electrochemical rate constant R^ d pro¬ 
duces approximately the same result as steady state voltammetry at a sphere on a surface of 
heterogeneous rate constant K® os = 7rln (2) A) ( ) ld /4, as was shown. 

We now investigate if this relation holds under the Marcus Hush formalism over a wide range 
of K° and A values. To do so, steady state voltammetry was simulated using the Marcus-Hush 
model described above for both a microdisc and a sphere on a surface. For each microdisc het- 
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erogeneous rate constant (K ^ d ) that was used, a corresponding sphere on a surface simulation 
was run with K® os = 7rln (2) K ° ld /4, with A the same in each case. Fig. [7] shows simulated 
steady state voltammograms at a microdisc electrode (solid line) and a sphere on a surface 
electrode (circles) at varying values of A (20, 24, 30, 40, 60, 100 and 200). K° is equal to 
5 x 10 -5 for the microdisc and (5 x 10~ 5 ) for the sphere. Qualitatively, the best agreement 
between the two geometries is seen at the smallest values of A where an almost exact fit is 
observed. Less good fits are seen at intermediate values of A. 

Fig. [8] (a) shows the percentage error between the fits in Fig. [3 as well as for three other 
values of K^ d : 5 x 10~ 3 (b), 5 x 10' 1 (c), and 5 x 10 1 (d). To obtain the corresponding values 
of AT s ° os , A)) ld is multiplied by pp. The percentage error is defined as: 

Error = Jsos ~ Jmd x 100% (54) 

J md 

It is seen that as K° increases, the error becomes less and less dependent on the values of A. 
This is because at high K° the Marcus-Hush model of electron transfer approaches the Butler- 
Volmer limit, and becomes independent of A. Over the whole range of K° and A values here 
studied, the error in the equivalence between a microdisc and a sphere on a surface never exceeds 
2.55% in regions of appreciable current. It is thus reasonable to say that the correspondence 
between the two geometries holds within the Marcus-Hush model of electron transfer, as well 
as within the Butler-Volmer model. 

We now apply our Butler-Volmer and Marcus-Hush models for voltammetry at a microdisc 
to the reduction of oxygen in two ionic liquids, and compare the effectiveness of the models in 
simulating experimental data. 
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4 Experimental 


4.1 Chemical reagents 

Ferrocene (Fe(C 5 H 5 ) 2 , Aldrich, 98%), tetra-n-butylammonium perchlorate (TBP, Fluka, Puriss 
electrochemical grade, 99%) and acetonitrile (MeCN, Fischer Scientific, HPLC grade, 99%) were 
used as received. Ionic liquids trihexyltetradecylphosphonium 6is(trifluoromethylsulfonyl)- 
imide ([P 14 , 6 , 6 , 6 ] [NTf 2 ]) and 1-butyl- 1-methylpyrrolidinium bis (trifluoromethylsulfonyl)imide 
([Bmpyrr][NTf 2 ]) were kindly donated by Professor C. Hardacre of Queen’s University, Belfast 
and were used as received. Oxygen (99.5%) was purchased from BOC, Surrey, UK. 

4.2 Instrumentation 

Electrochemical experiments (Cyclic Voltammetry (CV) and chronoamperometry) were con¬ 
ducted using a fi -Autolab potentiostat (Eco-Chemie, Netherlands). All experiments were con¬ 
ducted inside a temperature controlled Faraday cage.— The working microdisc electrode (either 
Gold or Platinum, IJ Cambria Scientific Ltd, UK), 10 /irn nominal diameter, was polished prior 
to use using a water-alumina slurry (1, 0.3, 0.05 /irn, five minutes on each grade) on soft lapping 
pads (Buehlcr, Illinois).— Precise radii were determined through calibration of the electrode 
with a 2.0 111 M solution of ferrocene in acetonitrile containing 0.1 M TBAP (silver wire as 
both a counter and quasi-reference electrode); chronoamperometry was recorded at 298 K, and 
assuming a diffusion coefficient of 2.3 x 10~ 9 m 2 s^ 1 ,— the data was analysed with respect 
to the Shoup and Szabo equation— (described in Section 12.711 . This gave electrode radii of 
5.10 ± 0.05 /irn (Gold) and 5.40 ± 0.05 /irn (Platinum). A 0.5 mm silver wire was used both 
as a counter and a quasi-reference electrode. The experiments conducted with a gold electrode 
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were arranged in a glass T-cell, as described previously,— with a plastic collar attached to the 
working electrode. A 10 /j,L aliquot of ionic liquid was transferred into the T-cell, which was 
left under vacuum (0.2 mbar) overnight, to ensure the complete removal of volatile substances, 
including atmospheric gases and any aqueous impurities. In the case of the platinum working 
electrode, a three electrode (with two silver wires acting as a counter and a reference electrode) 
system was used in a sealed sample vial with a pre-vacuumed ionic liquid. Oxygen was then 
introduced to the T-cell or sample vial via a drying column of molecular sieves (activated at 200 
°C overnight, 4 A, Sigma Aldrich). The gas was left on at a flow rate of 200 - 300 cm 3 min -1 for 
at least half an hour to reach saturation, at which point cyclic voltammograms were recorded 
to ensure equilibrium had been reached and the reductive current in successive experiments was 
stable. A wait time of minimum 15 minutes was imposed between each experiment to ensure a 
stable concentration of oxygen at the electrode surface. 

5 Results and discussion 

This section reports the comparative outcome of the use of BV theory and SMH theory to model 
experimental cyclic voltammetry of the reduction of oxygen in the ionic liquids [Bmpyrr] [NTf 2 ] 
and [P 14 , 6 ,6, 6 ][NTf 2 ]: 

0 2 + e ^ O* (55) 

This is followed by the inference of physical implications from the application of SMH and the 
determination of a value for solvent reorganisation energy, A, for both RTILs. 
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5.1 Reduction of Oxygen in [Bmpyrr] [NTfo] on a fi Au electrode 

The first system under study is the reduction of oxygen in [Bmpyrr] [NTf 2 ]. This ionic liquid 
was chosen for its relatively low viscosity (70 cP at 298 K)— and previously observed well 
characterised voltammetry for oxygen reduction.— Cyclic voltammetry (CV) was performed on 
saturated solutions of oxygen in [Bmpyrr] [NTf 2 ] over a range of scan rates (100-400 mV s _1 ). A 
blank CV, from -3 V - 2 V vs Ag, was recorded prior to the introduction of oxygen to ensure the 
ionic liquid was free from impurities, and to estimate solution capacitance for each scan rate. 
The capacitative current was subtracted from the experimental current to generate corrected 
voltammograms prior to further analysis. In all further experimental cases, the potential was 
swept from 0.00 V (vs Ag), a voltage at which ca. zero Faradaic current flows, to a potential 
cathodic of the reduction of oxygen (-1.50 V vs Ag), then back to 0.00 V (vs Ag). At a scan 
rate of 200 mV s^ 1 , the oxidative peak current at -0.8 V was 1.7 nA, the limiting current was 
-4.4 nA, and the half peak potential was -0.9 V vs Ag wire. The peak corresponding to the 
reduction of oxygen appeared on the forward scan at ca. -1.3 V vs Ag, whilst the oxidative 
peak appeared on the reverse scan at ca. -0.8 V vs Ag. The resultant voltammograms (Fig. 
El a-c) show more steady-state-like behaviour in the forward wave, which corresponds to the 
reduction of oxygen to superoxide, than in the back wave, corresponding to the reoxidation of 
superoxide to oxygen. The disparity in diffusion coefficients of the reactant and product in the 
ionic solvent results in different types of diffusion being observed at the electrode surface. In 
the case of oxygen, we see convergent Fickian diffusion of oxygen, whereas for the more slowly 
diffusing superoxide we observe planar Fickian diffusion.— 

CVs were modelled using a previously reported program— that utilises Butler-Volmer ki- 
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netics, as outlined above. As this is a well-reported system, 21 i 31 i 62 i 63 the starting point for the 


diffusion coefficient and the concentration of oxygen were taken from literature. Optimised 
values for k °, a, co 2 , Do 2 , D Q .-, and E'f (tabulated in Table [3]) were obtained to give the “best 
fit”, which was quantified with the calculation of a Mean Scaled Absolute Deviation (MSAD, as 
described in Section l2T8]l . These values were used across the scan rates with all values remaining 
constant, with the exception of Ef , which was allowed to deviate slightly (±0.02 V) to account 
for the use of a quasi-reference electrode. 61 The optimised values and the mean MSADs across 
the scan rates are reported in Table [3j 

The value of D 0 .-, the diffusion coefficient of the superoxide ion, is an order of magnitude 
smaller than that of the diffusion coefficient for oxygen, Dq 2 . This is due to the charged super¬ 
oxide ion’s increased interaction with the ionic solvent, which reduces its movement through 
the solution. This behaviour has been reported in other ionic liquids by Buzzeo et al ,— Evans 
et al,— and Huang et al— The value for the transfer coefficient, a, of 0.29, suggests an early, 
reactant-like transition state.- Another notable outcome is the value for k°, of 1.5 x 10~ 3 cm 
s -1 , which is comparable to k° for oxygen reduction in other ionic solvents, as seen in Tabled 
which details transfer coefficients and rate constants for oxygen reduction in a range of ionic 
liquids. 

Symmetric Marcus-Hush (SMH) theory was next used to model the electrode kinetics, and 
the same procedure as outlined above was followed to establish independent, optimal values for 
the parameters of concern. In this case, the transfer coefficient, a, is no longer a parameter 
that is taken into consideration. Instead, the reorganisation energy, A, is used. The standard 
electrochemical rate constant, k°, also takes on different values in this model. The best fit 
values and the mean MSAD values are outlined in Table [3J 
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The experimental value for A was determined to be 0.4 ± 0.1 eV in this solvent. This was 
determined by establishing best-fit parameters for the SMH model with experimental data. 
The BV program reproduces experimental data with a greater level of accuracy that the SMH 
program, as evidenced by lower MS AD values and a better fit, as seen in Table [3] and Fig. 
EJ This suggests BV kinetics parameterise the data better than the SMH model. A similar 
conclusion has been reached with respect to the modelling of redox processes in molecular 

solvents.JA^J, 

5.2 Reduction of Oxygen in [P 14 ,6,6,e] [NTf 2 ] on a /iPt electrode. 

The second system under consideration is the chemically irreversible reduction of oxygen in 
[P 14 , 6 , 6 , 6 ] [NTf 2 ]. This ionic liquid has high hydrophobicity— and previously reported well defined 
voltammetry.— Cyclic voltammetry (CV) was performed on saturated solutions of oxygen in 
[P 14 , 6 , 6 , 6 ][NTf 2 ] over a range of scan rates (100-400 mV s -1 ) (see Fig. (TUI a-c). A blank CV was 
recorded prior to the introduction of oxygen to ensure the ionic liquid was free from impurities. 
The baselines of the raw experimental CVs were adjusted to generate corrected voltammograms 
prior to further analysis. In all cases, the potential was swept from 0.00 V vs Ag, to a potential 
cathodic of the reduction of oxygen (-2.3 V vs Ag). At a scan rate of 200 mV s _1 , the limiting 
current was -6.5 nA, and the half peak potential was -1.6 V vs Ag. The reduction of oxygen 
occurred on the forward scan at -2.20 V vs Ag, whilst an oxidative feature appeared on the 
reverse scan at -0.86 V vs Ag. 

Unlike the case discussed in Section 15.11 in this situation the highly reactive superoxide that 
is generated by the initial reduction encounters a proton source in the form of the cation of 
the solvent. Evans et. al— detailed the mechanism of this process (shown in Scheme 1) and 
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Scheme 1 

suggested that the observed back peak in the experimental voltammetry could be an oxygen- 
derived species created in follow-up reactions, most likely HO*. Given the documented chemical 
irreversibility of this system, only the forward wave was considered for further analysis. Both 
the Butler-Volmcr and Marcus-Hush programs were modified to account for a homogeneous 
chemical step in the mass transport equations, as detailed in the Theory section. 

For the purposes of determining a starting point for Do 2 and concentration, potential step 
chronoamperometry was performed on this system by stepping from a potential at which no 
current flows, 0.0 V, to -2.0 V for 2s (see Fig. HU). This data was then analysed with respect to 
the Shoup and Szabo^ equation (as detailed in the Section 2.4), to simultaneously determine 
the concentration and diffusion coefficient of oxygen (see Table ED . 

Optimised values for k°, a , co 2 , Do 2 , and Ef were obtained to give the “best fit” using the 
modified Butler-Volmer program described above, which was quantified with the calculation of 
a Mean Scaled Absolute Deviation (MSAD, as described in the Section 12.811 . These values were 
used across the scan rates with all values remaining constant, with the exception of Ef, which 
was allowed to deviate slightly (± 0.02 V) to account for the use of a quasi-reference electrode. 
The optimised values and the mean MSADs across the scan rates are reported in Table [5j 
The values obtained for Do 2 and co 2 are consistent with previous literature reports and are 
in agreement with reported diffusion coefficients and saturated concentrations of oxygen in a 
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range of ionic liquids (Table |4j. 

The value for the transfer coefficient, a, of 0.31, suggests an early, reactant-like transition 
state. This is congruent with the result found in section 15.11 and a similar k° of 1 x 10~ 3 cm 
s~ x is reported, which is again in agreement with the literature values reported in Table SI 

SMH theory was used independently to generate the same set of parameters as generated 
in the BV theory with the transfer coefficient, a, once again being replaced by a representation 
for reorganisation energy, A. The experimental value for A was determined to be 0.5 ±0.1 eV 
in this solvent. 

Once again, the BV program provides parameters with a significantly greater level of cer¬ 
tainty that the SMH program. And, although the SMH theory gives greater physical insight 
compared with the phenomenological BV theory, this greater uncertainty leaves BV as the 
preferred theory, in terms of attaining more reliable parameterisation. 

5.3 Comparison of Oxygen Reduction in Pyrrolidinium and Phos- 
phonium based Ionic Liquids 

Both systems under study, oxygen reduction in either a pyrrolidinium or phosphonium based 
1L, display similarly low diffusion coefficients ( Do 2 = 2.05 and 3.95 xlO -10 m 2 s^ 1 respectively) 
and solubilities of oxygen (co 2 = 10 and 7.3 mM). In addition to this, both the transfer coeffi¬ 
cients (a = 0.29 ± 0.01, and 0.31 ± 0.01 respectively) and the reorganisation energies (A = 0.4 
± 0.1, and 0.5 ± 0.1 respectively) are similar. As outlined in the introduction, A encompasses 
contributions from both inner-sphere reorganisation of the molecule itself, and outer-sphere 
solvent reorganisation about the electroactive species. The inner-sphere reorganisation is due 
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to the change in bond length on going from O 2 to O* - (0.12 A),- with outer-sphere solvent 
reorganisation dominating in aqueous systems. Literature reports suggest the overall reorgani¬ 
sation energy for the oxygen self-exchange reaction in water is 1.26-4.51 eV.— The inner-sphere 
contribution was estimated to be 0.4-0.6 eV.- Given the close agreement with previously re¬ 
ported values for inner sphere reorganisation, 27, the experimental values for A in both of the 
systems considered in this work are consistent with a process dominated by an inner sphere 
contribution. Previous reports have indicated that oxygen itself had weak interactions with an 
ionic liquid,— and the result reported in the present work suggests that the ionic liquid does 
not undergo any significant reorganisation on moving from oxygen to superoxide. This may be 
related to the very small size of the C^/O* - species compared to the component ions of the 
RTIL. 

The reported best fit values for the rate constant given by SMH theory are similar to 
those established by BV theory, however there is less agreement with the experimental data 
than is seen when the BV program is used. This is reflected both in the larger mean MSAD 
values, and qualitatively observing the fit as seen in Figs. Eland [TO] Further studies may involve 
consideration of asymmetric MH theory, which relaxes the requirement that the potential energy 
parabola are of equal curvature, and introduces a new parameter to describe the disparity in 
curvature between the reactant and product. This has been considered in a number of molecular 
systems,and has been determined to fit experimental data more accurately than SMH. 
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6 Conclusions 


Models were developed to simulate cyclic voltammetry at a microdisc electrode and a sphere 
on a surface electrode using both Butler-Volmer and symmetric Marcus-Hush kinetics. The mi¬ 
crodisc model was used to investigate the possibility of kinetically limited steady state currents 
at a microdisc electrode, and it was found that this phenomenon is more likely to be observed 
at nano-sized sphere on a surface electrodes ( e.g . an impacting nanoparticle) because of the 
smaller size. It was also found that the steady state equivalence (subject to a mass transport 
correction) between microdisc electrodes and spherical electrodes resting on a surface also hold 
when Marcus-Hush kinetics are employed. 

The reduction of oxygen at a microdisc electrode in two room temperature ionic liquids, 
[Bmpyrr][NTf 2 ] and [P 14 ,6,6,6] [NTf 2 ], has been studied with cyclic voltammetry at 298 K. BV 
theory and SMH theory were used to create simulation programs, and the simulated voltam- 
mograms were compared with experimental data. BV theory, although empirical, provided 
voltammograms with lower MSAD values than those generated using SMH theory. A reorgan¬ 
isation energy, from the SMH program, was determined to be 0.4-0.5 eV for oxygen reduction 
in the two RTILs studied. This low value likely corresponds to a dominant inner sphere reor¬ 
ganisation, and limited solvent reorganisation. 

Further studies to investigate whether AMH theory is better than SMH in its ability to 
parameterise experimental data in ionic solvents will likely be valuable. 
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[ p i4,6,6,6l + [Bmpyrr] + 

Figure 1: RTILs used in this study; 

6is(trifluoromethylsulfonyl)imide [Pi 4 , 6 , 6 , 6 ][NTf 2 ] and 
6is(trifluoromethylsulfonyl)imide [Bmpyrr] [NTf 2 ]. 
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trihexyltetradecylphosphonium 
1-butyl-1-methylpyrrolidinium 
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(a) (b) 



Figure 2: (a) Schematic diagram of the microdisc electrode inlaid in a surrounding insulating 
material, defining the spatial coordinate system used in simulations, (b) The two dimensional 
simulation space used in this study. 
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Figure 3: Schematic diagram of the spatial grid used in microdisc simulations. 
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Figure 4: (a) Schematic diagram of an electroactive spherical nanoparticle resting on a conduc¬ 
tive, inert material, defining the spatial coordinate system used in simulations, (b) The two 
dimensional simulation space used in this study, (c) Schematic diagram of the spatial grid used 
in simulations. 
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Figure 5: Simulated steady state voltammetry microdisc electrode at various A values and K°. 
(a) K° = 5 x 1CT 1 (b) K° = 5 x 10~ 3 and (c) K° = 5 x 10~ 5 . In each case the dotted line 
represents the Butler-Volrner result (a = 0.5). A values are 20 (furthest from Butler-Volmer 
result), 25, 30, 40, 60, 100 and 200 (closest to Butler-Volmer result). 
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Figure 6: Schematic diagram of three electrode geometries which produce steady state voltam¬ 
metry described by Equation [50j (a) Isolated spherical electrode, (b) spherical electrode resting 
on a conductive, inert surface and (c) a microdisc electrode. 
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Figure 7: Comparison of steady state voltammetry at a microdisc electrode (solid lines) and a 
sphere on a surface electrode (circles) using Marcus-Hush kinetics at various values of A. For 
the microdisc electrode, K° = 5 x 10 -5 , and for the sphere K° = (5 x 10~ 5 ). A is equal to 

20 (bottom result), 25, 30, 40, 60, 100 and 200 
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Figure 8: Percentage error in the correspondence between steady state voltammetry at a mi¬ 
crodisc and a sphere on a surface for various values of K° and A. (a) Jl° ld = 5 x 10~ 5 , (b) 
Kmd = 5 x 10~ 3 , (c) K^ d = 5 x 10-\ (d) K^ d = 5 x 10 1 . Corresponding K® os values are 
obtained by multiplying Jl° ld by A values are indicated on the graphs. 
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Figure 9: Experimental cyclic voltammetry of the reduction of oxygen in [Bmpyrr][NTf 2 ] on 
a pAu electrode (solid line) at 298 K, compared to Butler-Volmcr Theory (dashed line), and 
symmetric Marcus-Hush Theory (circles) for scan rates of (a) 100 mV s—1, (b) 200 mV s -1 , 
and (c) 400 mV s~F Parameters corresponding to each simulation are detailed in Table [3j 
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Figure 10: Experimental cyclic voltammetry of the reduction of oxygen in [Pi 4 , 6 , 6 , 6 ] [NTf 2 ] on 
a pPt electrode (solid line) at 298 K, compared to Butler-Volmer Theory (dashed line), and 
symmetric Marcus-Hush Theory (circles) for scan rates of (a) 100 mV s -1 , (b) 200 mV s -1 , and 
(c) 400 mV s -1 . Parameters corresponding to each simulation are detailed in Table [5j 
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Figure 11: Representative experimental chronoamperometry of the reduction of oxygen in 
[P 14 , 6 , 6 , 6 ] [NTf 2 ] on a /iPt electrode (solid line) compared with fitting obtained using the Shoup- 
Szabo analysis (circles). 
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Tables 
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Parameter Description 

Units 

a 

Transfer coefficient 

Unit less 

Ci 

Concentration of species i 

mol m -3 

C? 

Bulk solution concentration of species i 

mol m -3 

C? 

Electrode surface concentration of species i 

mol m -3 

A 

Diffusion coefficient of species i 

2 —1 
ms 

E 

Applied potential 

V 

E f 

Formal potential of A/B couple 

V 

F 

Faraday constant = 96485 

C mol -1 

I 

Current 

A 

k 

(Pseudo) first order homogeneous rate constant 

s' 1 

k° 

Electrochemical rate constant 

m s _1 

k red/ox 

Reduction/oxidation rate constant 

Unit less 

A 

Reorganisation energy 

eV 

V 

Scan rate 

V s" 1 

R 

Gas constant = 8.314 

J K” 1 moU 1 

r 

Radial coordinate 

m 

r e 

Radius of electrode 

m 

T 

Temperature 

K 

t 

Time 

s 

z 

z coordinate 

m 


Table 1: List of symbols 
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Dimensionless Parameter Definition 


a 

D[ 

J 

K 

K° 

A 

R 

e 


O, 

D a 


2irFD A c\r e 


iLu 

D a K 


re UP 
D A K 


_F_ \ 
RT* 


RT 


(E - Ef) 


T 




Table 2: Dimensionless parameters. Species A refers to the species initially present in solution 
before the experiment /simulation begins. 
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Parameter 

Butler-Volmer 

Marcus-Hush 

Dq 2 / m 2 s" 1 

2.05 x 10~ 10 

2.05 x 10" 10 

D q .- / m 2 s^ 1 

1.80 x 10~ n 

1.55 x 10~ n 

co 2 / mM 

10 

10 

k° / cm s^ 1 

0.0023 

0.0015 

Ef /V 

-0.87 ±0.02 

-0.87 ±0.02 

a 

0.29 

- 

A / eV 

- 

0.4 

MSAD (average) / % 

9.52 

13.9 


Table 3: Simulation Results for the reduction of oxygen in [Bmpyrr] [NTf 2 ] on a /iAu electrode. 
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Ionic Liquid 

co 2 / mM 

Dq 2 / m 2 s" 1 

D or / nr> 

k° / cm s 1 

a 

[Emirn] [BF 4 ]— 

1.1 ±0.2 

1.7 ±0.2 x 10“ 9 

- 

0.94 ±0.13 x 10" 3 

0.46 ± 0.02 

[Pmim] [BF 4 ]^ 

0.97 ±0.05 

1.3 ±0.2 x 10" 9 

- 

1.5 ±0.4 x 10~ 3 

0.36 ±0.03 

[Bmirn] [BF 4 ]— 

1.1 ±0.1 

1.2 ±0.1 x 10" 9 

- 

0.83 ±0.2 x 10” 3 

0.36 ±0.02 

[MOPMPip] [NTf 2 ]— 

14.3 

1.0 x 10" 10 

- 

- 

0.18 

[HMPyrr] [NTf 2 ]— 

14.5 

2.5 x 10" 10 

- 

- 

0.13 

[Emirn] [NTf 2 ]— 

3.9 

8.3 x 10" 10 

1.27 x 10" 10 

- 

- 

[N6,2, 2 , 2 ][NTf 2 ]— 

11.6 

1.48 x 10~ 10 

4.66 x 10- 12 

- 

- 

[C 4 dmim] [NTf 2 ]— 

4.4 

4.7 x 10- 10 

7.5 x 10- 11 

0.94 ±0.13 x 10" 3 

0.4 

[N 6 , 2 , 2 , 2 ][NTf 2 ]— 

4.5 

4.1 x 10’ 10 

7.3 x IQ' 12 

0.3 x 10“ 3 

0.36 ±0.03 


Table 4: Literature values for Butler-Volmer fitting parameters for the reduction of oxygen in 
various ionic liquids 
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Parameter 

Butler-Volmer 

Marcus-Hush 

Do 2 / m 2 s _1 

3.95 x lO" 10 

3.95 x 10" 10 

co 2 / mM 

7.3 

7.3 

k° / cm s^ 1 

0.001 

0.001 

Ef /V 

-1.47 ±0.01 

-1.55 ±0.01 

a 

0.30 

- 

A / eV 

- 

0.5 

MSAD (average) / % 

2.80 

10.2 


Table 5: Simulation Results for the reduction of oxygen in [P 14 , 6 , 6 , 6 ] [NTf 2 ] on a yitPt electrode. 
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